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We present the first direct evidence for two-band superconductivity in the quasi-one dimensional 
chalcogenide Nb 2 Pdo. 8 iS 5 . Soft point contact spectroscopic measurements reveal Andreev reflection 
in the differential conductance G in the zero-resistance superconducting (SC) state below Tc(=Q.& K). 

Multiple peaks in G were clearly observed at 1.8 K and were successfully explained by the two- 
band Blonder-Tinkham-Klapwijk model with two gaps Ai = 0.48 meV and A 2 = 1.05 meV. Their 
evolution in temperature and magnetic field is consistent with the conventional BCS theory. The SC 
gap to Tc ratio for Ai and A 2 is 1.7 and 3.7, which is similar to the weak coupling BCS prediction 
of 3.5. These results demonstrate that the newly discovered niobium chalcogenide is a two-gap 
superconductor in the weak coupling limit. 


Materials with low-dimensional electronic structure 
have attracted interest because they provide a rich av¬ 
enue to explore various novel phases arising from strong 
electron correlations A spin (or charge) density 

wave is one such example that occurs due to geometry of 
the Fermi surface in momentum space that is conducive 
to the development of a Peierls instability which leads to 
a spatial modulation of the spins (or charges). Unconven¬ 
tional superconductivity with a triplet superconducting 
(SC) order parameter, or a chiral d-wave symmetry was 
also proposed due to strong repulsive electron-electron 
interactions An inhomogeneous superconducting 

phase with a periodic modulation in SC amplitude, or 
Fulde-Ferrel-Larkin-Ovchinnikov state, may also be re¬ 
alized at low temperature and high magnetic fields in 
low-dimensional compounds 0, [1| ■ 


The newly discovered chalcogenide superconductor 
Nb 2 Pda;S 5 (or Nb215) has attracted interest because of 
its quasi-one-dimensional (ID) electronic structure 


14| . Spin susceptibility is estimated to be strongly en¬ 


hanced at small q wavevectors, indicating that Nb215 is 
in close proximity to a magnetically ordered state with 
long modulation wave-lengths. When compared to its 
relatively low SC transition temperature (Tc = 6.6 K), 
its upper critical field is very high {Hc 2 = 37 T), exceed¬ 
ing the Pauli limiting field expected for weakly coupled 
superconductors {Hp 1.84Tc = 12 T). The tempera¬ 
ture dependence of the upper critical held perpendicular 
to the chain direction (b-axis) does not saturate at low 
temperatures, but linearly increases with decreasing tem¬ 
perature. When combined with the strong dependence 
on temperature of the Hc 2 anisotropy and multiple Fermi 
surfaces, the linear-in-T dependence of Hc 2 of Nb215 is 
indicative of multiband superconductivity Q. Recent 
muon spin relaxation and rotation measurements, how¬ 
ever, claimed a single s-wave energy gap from the tem¬ 
perature dependence of the magnetic penetration depth, 
thus requiring a more dehnitive study on the nature of 


SC gap in order to understand the SC mechanism of this 
low-dimensional superconductor. 

Here, we report the hrst, direct evidence for two su¬ 
perconducting gaps in the chalcogenide superconductor 
Nb 2 Pda;S 5 {x = 0.81). Differential conductance G{= 
dl/dV), which is obtained via soft point-contact spec¬ 
troscopy, shows an enhancement at zero-bias voltage be¬ 
low Tc due to the Andreev rehection and multiple peaks 
become more distinguishable with decreasing tempera¬ 
ture. A single band BTK (Blonder-Tinkham-Klapwijk) 
model fails to explain the spectroscopic features deep in¬ 
side the SC state, while the two band BTK model best 
explains the two peak structure with Ai = 0.48 meV and 
A 2 = 1.05 meV at 1.8 K. The magnetic held suppresses 
both the large and the small SC gaps, where a least- 
squares-ht with the BCS model reasonably explains the 
held dependence of both of them with Hc 2 = 17 T for 
H//c, which is consistent with the previous report [^. 
Taken together, these results prove the existence of mul¬ 
tiple superconducting gaps in the newly discovered quasi- 
ID chalcogenide compound Nb 2 Pd 2 ,S 5 . 

Single crystalline Nb 2 Pda;S 5 samples with x = 0.81 
were grown by solid reaction crystallizing in a monoclinic 
structure, where the Nb-S trigonal prisms are stacked 
along the b-axis forming conduct ing chain along that di¬ 
rection (see inset of Fig. 1) [l^[l^. Electrical resistiv¬ 
ity measurements were performed through the conven¬ 
tional four-probe technique on a needle like single crys¬ 
tal, where the electrical current is applied along the elon¬ 
gated b-axis. As shown in Fig. 1, the resistivity ratio 
between room temperature and T^, p{,(300 K)/pb{Tc), is 
approximately 3, which is similar to the one previously 
reported 0. The resistivity pb shows a sharp SC transi¬ 
tion to a zero value at 6.6 K. In the following, we analyze 
the data below 6.6 K, where the zero-resistance state is 
established. In order to perform the soft point-contact 
spectroscopy (SPCS), silver grains from Dupont 4929N 
were applied on the surface of the sample m- Magnetic 
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FIG. 1. Electrical resistance of Nb 2 Pdo. 8 iS 5 for the electric 
current applied along the b-axis. Inset : Crystal structure of 
Nb 2 Pdo. 8 iS 5 . Yellow, blue, and green symbols represent S, 
Pd, and Nb atoms, respectively. 


field is applied perpendicularly to the needle axis, i.e., 
H _L 6—axis. Tens of contacts, made on two different 
crystals from the same batch, showed consistent spec¬ 
troscopic signatures with respect to one another, thus 
ensuring reproducibility of the differential conductance 
data reported in this Letter. 


Figure [2][a) shows a representative plot of current / as 
a function of the bias voltage V (on the right ordinate) at 
1.8 K (solid line) and 9 K (dashed line). The nearly linear 
I — V characteristics for > A as well as the absence of 
dips in the differential conductance G indicate that the 
effective contact size is close to the Sharvin limit 17. M- 
Here A is the size of the SC gap. The dependence of G 
on the bias voltage of 1^ at 1.8 K is shown on the left 
ordinate of Fig. [SJa) (open circles). Multiple peaks in 
G, which are almost symmetric with respect to the zero- 
bias voltage, indicates two SC energy gaps for the Nb215 
compound. As shown in Fig. [5](b), the second derivative 
of the differential conductance spectra, cPG/dV^, reveals 
pronounced dips that correspond to the multiple peaks 
in G. Even though data are scattered, the first deriva¬ 
tive of the I — V curve at 1.8 K is consistent with the G 
obtained from an ac modulation technique (see Fig. SI, 
Supplementary Information), indicating that the multi¬ 
ple peak structure in G represents intrinsic SC properties 
of the Nb215 superconductor. 


In order to quantitatively understand the multiple SC 
gaps for Nb215, the multi-gap BTK model is used and the 
total conductance G = dl/dV is the sum of the multigap 
components with a weighting factor Wi for the ith compo¬ 
nent, or G = '^^uJiGi{V). The tunneling transparency 
of the point contact barrier Zi and quasi-particle life¬ 
time broadening effect F^ were parametrized in the BTK 
model for the conductance between normal and super- 





FIG. 2. (a) Soft point-contact spectroscopy of Nb 2 Pdo. 8 iSs. 
Differential conductance G at 1.8 K is shown as a function 
of the bias voltage on the left ordinate (open circles). The 
current-voltage (I-V) characteristic curves are shown on the 
right ordinate for T =1.8 K (solid line) and 9 K (dashed line), 
for which the sample is in the superconducting and in the 
normal state, respectively, (b) The second derivative of the 
differential conductance, A^G/AV^ at 1.8 K. (c) Normalized 
differential conductance of the Nb215 at 1.8 K. Gsym is the 
symmetrized differential conductance with respect to negative 
and positive bias voltages, i.e. Gsym = (G+ -|- G-)/{G{V = 
+5mV} + G(V = —5ml/)). Here G+ and G- is the dif¬ 
ferential conductance at positive and negative bias voltage, 
respectively. After symmetrization, Gsym was divided by the 
valued at -|-5 mV. Dashed and solid lines describe the least- 
squares fits from the single band and two band BTK models, 
respectively. 


conducting phases 17, 

§^iV, T) = Io^ [f{E - eV, T) - f(E, T)]a{E)dE 

/ OO -I 7^ rr 

^ ^^^^h\^^)a{E)dE. ( 1 ) 


Here, f{E,T) is the Fermi-Dirac distribution function at 
a finite temperature T and Iq is the prefactor that de- 
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FIG. 3. (a) Evolution of the soft point contact Andreev reflec¬ 
tion as a function of temperature, (b) Temperature depen¬ 
dence of the SC energy gaps obtained from d^G/dV^. BTK 
results at 1.8 K are added for comparison. Solid lines are 
estimated from the BCS theory with Tc of 6.6 K. 
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FIG. 4. (a) Distribution of SC energy gaps obtained from 
d^G/dV^. (b) Statistical counts of SC energy gap values for 
ten different contacts. Gaussian fit is used to determine the 
average SC gaps. 


pends on the density of states of the metallic tip and the 
SC material. The BTK integral kernel, cr{E), is written 
as 


a{E) = tn 


1 + Tat hiE)\'^ + {tn - 1) 

\i + iTN-ih{Er\^ 


( 2 ) 


The transparency of the barrier is given by tn = 
1/(1 -I- and the complex function "f{E) = {Nq{E) — 
l)/Np{E), where Nq{E,r) = {E + iV)/y/{E + iV)^ - A2 
and Np{E,r) = A/ y/{E + iV)'^ - 

Figure HKc) shows the normalized differential conduc¬ 
tance Gsym at 1.8 K, where G+ at positive voltage and 
G- at negative voltage is symmetrized and divided by 
the value at -1-5 mV, i.e., Gsym = {G+ + G_)/(G(= 
-|-5mV) -I- G(= —5mV)). A single band BTK model 
(dashed line) could explain the dip at zero bias voltage, 
but fails to explain the second dip at around ±0.7 mV. In 
contrast the two band BTK model (solid line) captures 
the characteristic features with Ai = 0.62 ± 0.018 meV, 
A 2 = 1.15 ± 0.012 meV and the weight of the first gap 
wi = 0.46. The amplitude of the two SC gaps is close 
to the value of the dip positions in the second deriva¬ 
tive of the differential conductance d^G/dV^, as shown 
in Fig. 2b. We note that a small F value of less than 
0.01 equally fits well the normalized conductance Gsym, 
indicating both a clean interface and a quasiparticle life¬ 
time that is close to infinite at low temperatures. The 
least-squares fit based on the two-band BTK model best 


describes the data with effective tunneling barriers Z of 
0.3 and 0.5 for the small and large SC gaps, respectively. 

Figure EJa) selectively shows the differential conduc¬ 
tance G of Nb215 for several temperatures in its SC 
state. The conductance dip at the zero-bias voltage, 
which is typical for Andreev reflection when Z 7 ^ 0, 
becomes shallower with increasing temperature and dis¬ 
appears completely for temperatures above 2.9 K due 
to thermal smearing effects. Another dip structure at 
±0.7 mV for 1.8 K, which reflects the two-gap nature 
of the dichalcogenide superconductor, is also smeared 
out with increasing temperature due to thermal effects. 
The decrease in the dip position with increasing temper¬ 
ature is consistent with the decrease in the gap ampli¬ 
tude. The temperature dependence of the two SC gaps 
is shown in Fig. EKb), where the gap amplitudes are ob¬ 
tained from the two-band BTK fitting (solid symbols) 
and from the dip in the second derivative of the differ¬ 
ential conductance (open symbols) for several different 
contacts made on the ac-plane of different crystals. Solid 
lines depict the dependence expected from the weak cou¬ 
pling BCS theory, which reasonably describes the tem¬ 
perature dependence of the gaps with Ai = 0.48 meV 
and A 2 = 1.05 meV. 

Figure |4] shows statistical counts for ten different point 
contacts and the distribution of SC energy gaps obtained 
at 1.8 K from the second derivative of the differential con¬ 
ductance, d^G/dV^, for Nb215. The small SC gap Ai is 
concentrated between 0.4 and 0.6 meV, while the large 
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FIG. 5. (a) Magnetic field dependence of the differential con¬ 
ductance of Nb215 at 1.8 K for 0, 1, 3, 6, and 9 Tesla. The 
conductance for each magnetic field is vertically displaced for 
clarity, (b) Magnetic field dependence of the SC gap from 
d^G/dV^. Dashed lines correspond to least-squares fits from 
Eq (3). 


SC gap A 2 is between 1.0 and 1.15 meV. As shown in 
Fig. 4(b), the statistical distribution can be explained by 
a Gaussian function, where the peaks are located at 0.48 
±0.02 and 1.05 ±0.04 meV for Ai and A 2 , respectively. 
The SC gap ratios against Tc, 2Ao/kBTc, are 1.7 and 3.7, 
where the gap ratio for the large gap is close to the BCS 
prediction of 3.5 for weak coupling superconductors [21 1. 
Recently, the temperature dependence of the specific heat 
for Nb 2 PdS 5 was analyzed in terms of multiple SC gaps, 
where the small and large gap ratio is 1.9 and 6.4, respec¬ 
tively 22|. The difference between this work and Ref 221 
may arise from the limited temperature range over which 
the specific heat was studied, i.e., 2 K< T < 6.6 K. When 
measured down to 0.5 K, a similar specific heat study for 
another chalcogenide Nb 2 Pdi. 2 Se 5 showed that the gap 
ratios are 2.1 and 4.1, which is similar to this work [l6|. 
We note that the gap values obtained from the two-band 
BTK model are consistent with those from . 

The dependence of the differential conductance of 
Nb215 at 1.8 K on the magnetic field H is shown in 
Fig.[S)a), where H is applied perpendicularly to the crys¬ 
tallographic 6-axis. For clarity, the conductance for each 
field is vertically displaced with respect to the zero-field 
spectrum. The multiple peaks observed at zero field are 
rapidly suppressed with increasing field, which could be 
due to a fast decrease in the Andreev reflection from the 
small SC gap Ai under magnetic field. The dip observed 


at zero-bias voltage, which appears in the Andreev re¬ 
flection for a finite Z, becomes blurred with increasing 
magnetic field. The pair-breaking effects induced by the 
magnetic field are indicated by the enhancement in the 
effective broadening parameter P due to a reduction in 
the quasiparticle lifetime ItI 23|. 

The magnetic field dependence of the SC energy gaps 
as determined from the dip positions in d^G/dV"^ is 
shown in Fig. [5)b). The dashed lines are least-square- 
fits of the SC gaps for type-II superconductors, where 
the SC gap ratio 2AjknTr is assumed to be independent 
of magnetic field [2J, [25 1: 


A(r,R) = A(T,0)x/l-S/Rc2(0), (3) 

where A(r,0) = A(0, 0)tanh(l.78 VPc(S)/T - 1) and 
Tc{B) = 'Fc(0)\/1 — B/Bc2{0)- The upper critical field 
Bc 2 at 1.8 K is 17 T, which is consistent with magneto¬ 
resistance measurements [^. 

To conclude, we have directly shown that the quasi-lD 
chalcogenide Nb 2 Pdo. 8 iS 5 is a two-gap superconductor, 
where the amplitude is 1.05 and 0.48 meV for the large 
and the small SC gaps, respectively. The differential con¬ 
ductance curves measured from soft point-contact spec¬ 
troscopy for Nb215 were successfully explained by the 
two-band BTK model with a finite point contact barrier 
Z and a quasiparticle lifetime broadening effect F. With 
increasing temperature, the multiple peak structure and 
the dip at the zero-bias voltage in the differential conduc¬ 
tance become weaker due to thermal broadening. The 
temperature dependence as well as the magnetic field de¬ 
pendence of the SC gaps are in agreement with the BCS 
theory. The gap/temperature ratio for the larger gap is 
3.7, which is slightly larger than BCS value of 3.5 for 
weakly coupled superconductors. Although the isotropic 
BTK model was successfully applied to describe the dif¬ 
ferential conductance of the Nb215 superconductor, its 
quasi-lD Fermi surfaces have yet to be taken into ac¬ 
count to make a better understanding of the point con¬ 
tact spectra. 
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